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GAMMA RAY ASTROPHYSICS* 

F. W. Stecker 
Theoretical Studies Group 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 20771 

1. INTRODUCTION 

In the context of an advanced study institute on the origin of cosmic rays, 
one can look at the subject of theoretical gamma-ray astrophysics from two points 
of view. The first seeks to answer the question, "What is the origin of the ob- 
served cosmic gamma radiation? ", considering this radiation as a component of 
cosmic radiation. The second seeks to answer the question, "What does the 
observed cosmic gamma radiation tell us about the origin of cosmic rays?" 
with the term "cosmic rays" meant to be the primary nuclear component of 
cosmic radiation (or perhaps the electronic component) and the gamma radiation 
considered to be a secondary product of various interactions between primary 
radiation and fields, photons and nuclei in the cosmos. This latter question, in 
which most, of you here I'm sure are "primarily" interested, is the historically 
older question, but it is inherently a much more difficult one to answer and is 
one which cannot be unambiguously answered at this point in time, given our 
relatively primitive observations and the context of our present theoretical 

Lectures presented at NATO Advanced Study Institute on“The Origin of Cosmic Rays" Aug. 26- 
Sep. 6, 1974. 
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understanding. Indeed it may well be that the gamma ray background radiation 
is not even a product of cosmic ray interactions, at least up the energies thus 
far observed, and we will discuss this possibility later. Thus, of necessity, we 
will spend most of our time here discussing the origin of the gamma radiation 
itself, bearing in mind the relation of this problem with that of the origin of the 
nuclear component of cosmic radiation. 

The possible existence of a secondary component of cosmic gamma radiation 
first appeared in the literature almost incidental to problems bearing on the 
question of the origin and propagation of the primary cosmic radiation. After it 
had been determined that overwhelming component of primary cosmic radiation 
impinging on the upper atmosphere consists of high-energy protons, Feenberg 
and Primakoff (1947) addressed themselves to an explanation of the conspicuous 
absense of electrons in significant quantities in the primary radiation. To do 
this, they examined the various interactions which cosmic ray electrons and 
protons could be expected to undergo with low-energy starlight photons in inter- 
stellar space and found that electrons could be effectively depleted of their 
energy by Compton interactions in a fraction of the age of the universe. In this 
process, the energy of the electrons is transferred to the photons which can be 
boosted to X-ray or gamma ray energy. Similarly, Hayakawa (1952), in examining 
the propagation of cosmic radiation through interstellar space, pointed out the 
effect of meson-producing nuclear interactions between cosmic rays and inter- 
stellar gas. Hayakawa noted that the neutral poins produced would decay to 
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produce cosmic gamma radiation. The production of bremsstrahlung radiation 
by cosmic rays was discussed by Hutchinson (1952). However, the idea of 
establishing a science of gamma ray astronomy itself as a tool for answering 
questions in high-energy astrophysics and cosmology appears to have been 
stimulated in an important article by Morrison (1958). Morrison (1958) and 
Felten and Morrison (1963) pointed out that the processes of most significance 
for producing cosmic gamma rays were (1) electron bremsstrahlung, (2) electron- 
photon interactions (Compton effect), (3) cosmic-ray produced pi on decay, (4) 
synchrotron radiation, (5) annihilation produced pion decay, and (6) line emission 
from electron-position annihilation and nuclear deexcitation. 

We will first review here these various processes for producing cosmic 
gamma radiation, as well as the significant processes for absorption of cosmic 
gamma radiation, discussing the basic physics of these processes. We will then 
attempt to place these processes in their astrophysical context in the galaxy and 
the universe as a whole. We will then turn to the interpretation of the present data 
on cosmic gamma radiation and its implications for cosmology and cosmic ray 
origin. 

2. GAMMA RAY PRODUCTION PROCESSES 

A. Compton Scattering Between Cosmic-Ray Electrons and Starlight 
and Microwave Background Photons. 

Compton scattering is the relativistic limit of the electron -photon scattering 
process which is referred to as Thomson scattering in the non- relativistic case. 


3 



The traditional Compton effect is one in which a 7 -ray scatters off an electron 
at rest in the laboratory. In the case of astrophysical interest, the electron is a 
cosmic-ray of energy E = ymc 2 and velocity v = /3c. A photon of initial energy e 
then scatters off the electron, coming out of the interaction with an energy e'. 
Denoting quantities in the electron rest system by asterisks and denoting the 
scattering angle of the photon by u*, the initial and final energy of the photon in 
the electron rest system are related by 



1 + — — (1 - COS V*) 
me 2 


( 1 ) 


Denoting the angle between the electron and the photon by a, it then follows 

that 


e* = ye (1 + /i cos a) 


e' = ye*' (1-/3 cos a*') 

and 


so that 


tan a 


s in a 

7 (cos a + /3) 


( 2 ) 


y 2 e(l + /3 cos a) (1 - j3 cos a*') 

1 + (1 + fi cos a) (1 - cos v*) 

me 2 


(3) 
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and, on the average, the final photon energy in the observers frame is 


e' ~y 2 e (ye << me 2 ) (4) 

For the astrophysical conditions we will primarily be concerned with here, the 
condition y e < < me 2 holds and the cross section for the scattering is 


o- -* cr = ^.77 j J: | = 6.65 x 10 25 cm 2 

c T 3 \ 2 j 

\mc * I 


(5) 


(Heitler 1960). 

In the other extreme, the Klein-Nishina form of the cross section holds, 
which has the asymptotic form. 


°c “* 7T 


'eM 2 /mcA 

ilBC 2 / ' £ ' 


UlnfHl 

\ \mc J 


(ye > > me 2 ) 


( 6 ) 


Also in this extreme 


€ > ^ -y mc 2 ~ E 

However in most cases, i.e. whenever e E y « 
y -ray produced is on the average 


(7) 


(me 2 ) 2 , the energy of the 



( 8 ) 


In the important case of astrophysical interest where the cosmic ray electrons 
are assumed to have an energy distribution of power-law form 


I e (E) = KE" r 


(9) 
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if we make the delta-function approximation for the differential production 


function 


a(E y |€, E) *^8 (ly-l <e> 7 2 


(10) 


the y -ray production rate as a function of energy is given by 


q(V = 47Tn ph o- T /dEKE _r 5 (* y - i <e> 7 2 ) 



8 77 

“ "3~ °T^ph 


(mc 2 ) 1_r 





where n ph is the number density of target photons of low energy in the inter- 
stellar or intergalactic medium (usually in the form of starlight photons of ~ 1 eV 
energy or universal microwave blackbody photons of energy ~ 0.6 x 10~ 3 eV) and 
p ph = n ph <e> is the energy density of the radiation. 

In particular, for blackbody radiation of temperature T 



3.6 kT - 3. 1 x 10" 10 T MeV 
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and 


( 12 ) 


/ o ph = 4.7 5 x lCf 9 T 4 MeV/cm 3 

B. Synchrotron Radiation 
Synchrotron radiation, or magnetic bremsstrahlung, is the radiation emitted 
by a relativistic particle spiraling in a magnetic field. Its mathematical descrip- 
tion has been given by Schwinger (1949). An electron suffers energy losses by 
synchrotron radiation at a rate 

(H) = 

synch 

where the magnetic energy density p H is equal to H 2 /8 it. This rate can be com- 
pared with the energy loss suffered by an electron through Compton interactions 
with a photon field of energy density p ph . That rate is 

1 ^) = - ia cp (14) 

\dt / 3 T 

C 

The equivalence of equations (13) and (14) can be shown to be a direct consequence 
of electromagnetic theory (Jones 1965). The photons emitted as synchrotron 
radiation have a characteristic frequency given by 



The synchrotron effect can be thought of as the interactions of an electron 
with "virtual" photons of the magnetic field. Denoting these photons by the 
subscript v, it is found that they have an average energy 
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8 

Stt/3 


and a number density 


v 877 £ 


{Jones 1967). 


The resultant 7 -rays have a characteristic energy found from equation (15) 


to be 


— y 2 (±) 
2 \mc/ 


Because this is the same type of energy dependence as that for Compton 
radiation, a power-law cosmie-ray electron spectrum of the form KE will 
again generate a synchrotron radiation spectrum of the form 


q(V wE r 


The above equations can be used to find some useful numerical relations 
involving 7 -ray production by synchrotron radiating electrons, with the syn- 
chrotron radiation being determined from radio observations (e.g., in the 
10-1000 MHz range) and the 7 -radiation being produced from Compton inter- 
actions of the same electrons responsible for the synchrotron radio emission. 
Typical galactic magnetic fields are of the order of a few ^G. From equation 
(15) it follows that the frequency of a synchrotron radiating electron v as a 
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function of electron energy is given by 

v -4.2 y 2 H, MHz (20) 

with the magnetic field strength given in gauss. It then follows from equations 
(12) and (20) that the energy of a Compton y-ray produced by an electron which 
synchrotron radiates at frequency v s is given by 


E y = | CO 7 2 = 0.74 x 10' 10 v s MeV ' (21) 

To determine the relative importance of synchrotron radiation and Compton 
scattering in producing y -radiation under astrophysical conditions, we note 
that the energy of a synchrotron photon is 

E = hv = 2.8 x 10" 15 y 2 H, MeV (22) 

y % s s J* 

We then specify that an electron of energy y c me 2 will radiate a Compton photon 
at the same mean energy as an electron of energy y s me 2 will radiate a syn- 
chrotron photon. It then follows from equations (8), (12) and (22) that the ratio 



(23) 


The relative production rates from synchrotron and Compton radiation are then 
related by 

Q s _ ^h_ (24) 

" £ P h W 
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In the galaxy, for example, Pph - J , but y s » y c so that \ (yj « 1* ( y ) 
and synchrotron radiation is negligible compared to Compton scattering as a 
7 -ray production mechanism. 

C. Bremsstrahlung Interactions 

Bremsstrahlung, which is a German word meaning "braking radiation" is 
the radiation emitted by a charged particle accompanying deceleration. The 
cross sections for 7 -ray production from bremsstrahlung are derived in Heitler’s 
book (Heitler 1954). In the case of bremsstrahlung from non-relativistic elec- 
trons radiating in the field of a target nucleus, the cross sections for production 
of a 7 -ray of energy E y by an electron of energy E in the field of a nucleus 01 
charge Z is given by 


^(EjE) »a b (E, E y ) f(E r |E) . < ! 

where f b (E y |E) is the normalized distribution function for y -ray production, 
and where, in terms of the kinetic energy T = E - me 2 


°b (E> E r ) = 



(=£) m 



(26) 


for 2ttZcl « j3 « 1 


The cross section for bremsstrahlung by relativistic electrons is given by 


(E, E r ) = 


4a.Z 2 

77 


cr In 



— 1 ^ for me 2 « E « a" 1 Z‘ 1/3 mc 2 


(27) 
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In the ultrarelavistic case, where the cross section is calculated by taking 
into account the screening of the charge of the atomic nucleus by atomic electrons, 
the resultant cross section is given by 


ME. v = ^ ( 183 2-1/3 ^)[r r 


( 28 ) 


for E » a " 1 Z“ 1/3 mc 2 

Of course, in the case of an ionized gas, equation (27), which holds for the 
case when screening effects are unimportant, is applicable at all relativistic 
energies. The cross sections given in equations (26) and (27) may be corrected 
for additional contributions from interactions between cosmic ray electrons and 
atomic electrons by the replacement Z 2 - Z(Z + 1). 

Another case which may be of astrophysical interest for gamma-ray pro- 
duction may be that of bremsstrahlung from nonrelativistic protons, as has been 
pointed out by Boldt and Serlimitsos (1969). In this case, it is the electrons 
which is at rest and the proton which has the kinetic energy in the laboratory 
system (observers system). The appropriate cross section for this interaction, 
which corresponds to equation (26) for cosmic ray electrons, is then 


- b <V V = ^°-r 


M c 2 ^ 

_p 

r 

pr 


In- 


T pr + v / (M g /M p ) T pr ~-TE;[/(M e /M p )T p ^ (29) 


E„ 




It is immediately obvious from equations (27) and (28) that for relativistic 
particles the bremsstrahlung cross sections have little or no dependence, except 
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for a linear one, on E. Indeed, equation (28) may be written in the form 


^, E r >--§(i) m 

where <(M)>is the average mass of the target atoms in grams and 0^> is the 
average radiation length for the gas in grams per cm 2 . The average radiation 
length for interstellar matter is 


X = 65 g/cm 2 ' 01 

based on the values given for pure hydrogen and pure helium by Dovshenko and 
Pomanskii (1964) of 


Xjj = 62,8 g/cm 2 


(32) 


and 


^ = 93.1 g/cm 2 


(33) 


To a good approximation, especially in the case of relativistic bremsstrahlung 
(see Heitler 1954) the normalized distribution of gamma-rays produced may be 
taken to be a square distribution given by 

fE" 1 for 0 < E y < E 


f (E r |E) = < 


(34) 


\0 otherwise 

so that the gamma-ray production spectrum is given by 
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^y) = 


<M> 

<X> 


it 

L E r 


dEI (E) 




( 35 ) 


<X> J ' E 


where p( r) is the matter density of the gas in grams per cm . 
between cosmic-ray electrons and interstellar gas we may use 
write equation (35) in the form 


For bremsstrahlung 
equation (31) to 


q. (E ) = 4.3 x 10~ 2S n Ie( ' Ey ' ) cm' 3 s ' 1 MeV" 1 (36 > 

y 

where n is the number density of nuclei in the production region. 

D. Cosmic-Ray Produced tt° Meson Decay 

f 

Cosmic-rays of high enough energy can produce various secondary particles 
of short lifetime upon collision with interstellar or intergalactic gas nuclei. Of 


these secondaries, the most important for production of y -radiation is the tt° 
meson which decays almost 100% of the time into 2 gamma rays. An extensive 
treatment of these various secondary production processes leading to cosmic 
y-ray production has been given elsewhere (Stecker 1970, 1971a, 1973a) and the 
reader is referred to these references for more detail of the material outlined 
here. 

The collision processes of highest frequency occurring in interstellar and 
intergalactic space are those between cosmic-ray protons and hydrogen gas 
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nuclei and are therefore p-p interactions. The threshold kinetic energy which a 
cosmic -ray proton must have to produce a secondary particle of mass m in such 
an interaction is 


T 


th 



m 


(37) 


\ p/ 

Thus, to produce a single tt° meson of rest mass 135 MeV/c 2 requires a thresh- 
hold kinetic energy of 280 MeV/c 2 . Within the past two decades, many measure- 


ments have been made of the cross sections for inclusive production of various 


secondary particles in high energy interactions using proton accelerators. These 
data are summarized in Figure 1, taken from Stecker (1973a) wherein the references 
may be found. Figure 1 shows the product of the production cross section 0 
and multiplicity C n 0 given as a function of kinetic energy T. Utilizing the data 
shown in Figure 1, and a demodulated cosmic-ray proton spectrum given by 
Comstock, et al. (1972), Stecker (1973a) calculated a totaly -ray production rate 
from this process in the solar region of the galaxy to be 


q y - (1-3 ± 0.2) x 10“ 2S r» H cm -3 s -1 (38) 

p 

taking account of the effect of p-He and a -H interactions as well as that of p-p 
interactions. 

Figure 2, which shows the product of the cosmic-ray energy distribution 
I(T) and the function cr Q £ Q of Figure 1, is a measure of the kinetic energy of 
the cosmic -ray protons which are most effective in producing w 0 mesons in the 
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T(GeV) 


Figure 1. Cross section times multiplicity for neutral pion production in p-p interactions 
as a function of incident kinetic energy (Stecker 1973a). 








galaxy. It is found from this figure that the average proton which produces a rr° 
meson in the galaxy has an energy of about 3.3 GeV and that 85% of the 77 0 mesons 
produced in the galaxy are produced by cosmic-ray protons in the energy range 
between 1 and 10 GeV. 

When a 77 ° meson decays into 2 7 -rays, because of conservation of momen- 
tum they both are produced with an energy of m^ c 2 /2 = 67.5 MeV in the rest 
system of the pion. However, in the system of a terrestrial observer, in general, 
these 7 -rays have unequal energies which are determined by the relativistic 
Lorentz transformation. To determine the observed 7 -ray energy, let 7 ^ and/3^ 
refer to the energy and velocity of the parent 77 0 meson in the observer's system. 
Let us also define the energy 1 / = c 2 /2. Then the energies of the 7 -rays in 

the observer's system are given by the Doppler relation 

E r i,2 = ^C 1 ± A* cos e ) (39 ^ 

where 0 is the angle between the 7 -ray and the axis of the transformation in the 
pion rest system. The ± sign applies because the 7 -rays come off at opposite 
directions in the pion rest system. 

Since there is no preferred direction to the decay in the pion rest system, the 
distribution in cos 8 is isotropic, i.e., the distribution function 

f(cos O') = const. 

and therefore , for a given value of y* A ' © 

f(E r |E„) = const. (41) 
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The extreme upper and lower limits on the 7 -ray energy in the observers system, 
E u and respectively, are defined by letting cos 6 = 1 in equation (39). We 
then find 


E u = V yA l + Ar) 

and 

ffC^^d-Ar) 

and since = (1 - /3 2 ) 1/2 , it follows that E u = v 2 or 

In E + In Ejj 

1 u 4, 

In v - 

2 


(42) 


(43) 


From equations (41) and (42), it follows that the normalized 7 -ray energy 
distribution function from the decay of pions of energy E is given by 


(E^-M 2 c*) 


-1/2 




for Ec < E_ < E 
4- y u 


otherwi se 


(44) 


and thus the total production rate from cosmic -ray protons with an energy 
spectrum I(E p ) is 


I 


q(E y ) = 47 m dE p I (E p ) ■ 2£(E p ) 


“th 


(45) 


"Err 

", max 

x 1 dE^ 


f 
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To find the limits on the E integration in equation (45), we note that the 
relation between the maximum y- ray energy and E n can be written as 


E = E + __ 

77 U + E 


( 46 ) 


and we may therefore reverse the criteria to note that y -rays of energy E y 
may be produced by pions of energy as low as 


m in " ^7 + -g- 


(47) 


In the upper limit, we note that as - oo , -» 1 and from equation (42), 

Ep - 0. This implies that as we increase the pion energy without limit, the 
allowed range of y -rays produced expands in the lower limit to allow for all 
energies and pions of arbitrarily high energy can produce y-rays of energy E y , 
i.e., 


E - c° 

77, max 

and therefore equation (45) becomes 


(48) 


q(E y ) = 87m 


JE th 


dE p I (E p ) C(E P ) 


I 


dE, 


r o . 

E y +vVE y 




(49) 


It follows from equations (43) and (44) that the y -ray spectrum from an 
arbitrary distribution of energies in a decaying pion beam can be looked upon 
as the superposition of a number of rectangular distributions centered about ‘in v 
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on a logarithmic plot. Also, because the energy v is included within the allowed 
range of y - ray energies, there is a maximum in the energy spectral distribution 
at this value. Intuitively, one can also see from equation (43) that the full spec- 
trum from any arbitrary pion energy distribution should be symmetric about the 
maximum at ■In v on a logarithmic plot. These characteristics can be rigorously 
proven (Stecker 1971a). The differential y - ray energy spectrum calculated for 
galactic cosmic-ray interactions per unit path length and target nucleon density 
is shown in Figure 3 calculated from a two-component kinematical model for pion 
production in high-energy cosmic ray interactions (Stecker 1970). The integral 
y -ray energy spectrum is shown in Figure 4 along with calculations by Cavallo 
and Gould (1971) based on a different kinematical model. These spectra have been 
normalized to compare the shapes obtained. The wiggles in both spectral calcu- 
lations represent artifacts of the kinematical models assumed and should not be 
taken too seriously. The shapes of the two spectra are in good agreement and 
probably represent an accurate approximation to reality within the uncertainty 
indicated by the wiggles. 

Various other secondary products of high-energy cosmic-ray interactions 
such as other mesons, hyperons and excited nucleon states contribute to the 
overall cosmic y - ray spectrum, particularly at high energies. However, their 
total contribution to the production rate is relatively minor. The kinematics of 
these reactions is, in general, much more complicated than in the case of pion 
decay and has been thoroughly discussed elsewhere (Stecker 1971a). Figure 5, 
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i-COMPONENT (p-p) 
f-COMPONENT (p-p) 
TOTAL p-p 
TOTAL 

p-p+p-a+a-p+a-a 


Figure 3. The calculated differential production spectrum of y-rays produced in cosmic-ray 
interactions in the galaxy based on the "isobar <i) -plus -fireball (f) M model of Stecker {1970). 
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Figure 4. A comparison of the shapes of the integral galactic pion-decay energy spectra 
calculated by Stecker (1970) and Cavallo and Gould (1971). The total production rate is 
normalized to unity. 
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Figure 5. Total calculated galactic y-ray production spectrum 
from cosmic-ray interactions (Stecker 1971a). 





taken from Stecker (1971a) shows the total y-ray spectrum calculated for all 
secondary particles produced by interactions of cosmic rays with interstellar 
gas in the galaxy. 

E. Nucleon -Antinucleon Annihilations 

If antimatter exists in significant quantities in the universe, secondary 
particles can result from the annihilation of nucleons and antinucleons. Of these 
secondaries, the most significant for y-ray production again are the tt ° mesons. 
Thus, the resultant y -ray production spectrum is also symmetric about 67.5 MeV 
on a logarithmic plot, but in this case it is also bounded by the limited rest-mass 
energy which can be released in the annihilation. Frye and Smith (1966) using 
accelerator data, and independently Stecker (1967, 1971a) using a theoretical pion 
production model for nucleon-antinucleon annihilation have calculated the resultant 
y-ray spectrum from annihilations at rest. There is excellent agreement between 
the two calculations and the resultant spectrum is shown in Figure 6. 

F. Form of the Spectrum from Pion Production at High Energies 

It is interesting to examine the asymptotic form of the y-ray spectrum to 
be expected from the decay of 77 0 mesons produced at energies above a few GeV 
by cosmic rays having a power-law spectrum <* E”^. In this case, we assume 
that the pion production cross section is constant and that the pion multiplicity 
rises as a power of the primary energy <* E a . We also assume that the average 
pion energy rises as . We can then write the production function in the form 

^(EjE p )=a- 0 E«S(E 7r -X 0 Ep (50) 
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where the coefficients o 0 and x 0 and the exponents a and b are taken to be 
constants. Equation (49) then reduces to 


q(Ey) = 877nK p cr Q 


= 877nK P °o 




5 ( E 7T ~ Xo E p) 
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dE p 


f "[(F +b)-a] 



(51) 
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K cr„y Ug /b )“l] E“g /b 
p 0*0 y 


where 

g s [(r + b) - (a + 1)] 

As can be seen from equation (51), the high-energy y -ray spectrum produced 
by cosmic rays having an energy spectrum with index T is also a power-law, but 
one which is different and, in general, steeper than the primary spectrum. For 
example, in the case where F = 2.5 and a = 1/4, b = 3/4 (the Fermi model for 
meson production), they -ray index g/b = 2.67. The primary and y-ray spectra 

■7 

have the same index in the important case where a constant and large fraction 
of the primary energy of the cosmic ray is carried off by a single secondary 
which decays to produce a 77 ° meson. This corresponds to the case a = 0, b = 1 ■ 
and therefore g/b = T . 
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G. y -Rays from Photomeson Production at Ultrahigh Energies 
A process which may be important in producing y-rays of energies in the 
10 19 eV range is that of photomeson production, i.e., reactions of the type 

y + p - tt° + p 

It is now generally believed that the universe is filled with thermal micro- 
wave radiation of temperature 2.7 K which is a remnant of the primordial big- 
bang. This temperature corresponds to an average photon energy of 2.7kT = 

6.4 x 10 -4 eV and an average photon number density of about 400 photons per cm 3 . 
To an ultrahigh energy cosmic ray of energy in the 10 2 ° eV range, these micro- 
wave photons look like y -rays of energy equal to hundreds of MeV. These 
energies are above the threshhold for reaction (52) so that various ultrahigh 
energy mesons can be produced. One very important effect of this is the attenu- 
ation of ultrahigh energy cosmic rays in inter galactic space in 10 8 -10 9 years 
(Greisen 1966, Zatsepin and Kuz 'min 1966, Stecker 1968a). A typical y -ray 
produced by reaction (52) carries off about 107. of the primary energy of the 
cosmic ray (Stecker 1973b). The y-rays produced may themselves be attenuated 
by pair production processes of the type 

y 2-7 + y - e+ + e " 

and this process, together with Compton scattering of the electrons and positrons 
by the 2.7 K radiation may lead to a cascade process which has been treated in 
detail by Wdowczyk, et al. (1972) and Stecker (1973b). 
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3. y ~ Ray Absorption Mechanisms 

Various processes are of astrophysical importance in depleting 7 -rays in . 
the galaxy and the universe. By ’'absorption" we will mean not only those proc- 
esses in which the 7 -ray completely disappears, such as process (53), but also 
those processes in which the 7-ray is scattered out of the energy range of interest 
as can occur in the case of Compton scattering. 

We will consider two basic catagories of absorption processes: absorption 
in matter and absorption through interactions with radiation. The later process 
is of importance because of the existence of the 2.7 K thermal universal radiation 
field. 

A. Absorption through Interactions with Radiation 

Let us first consider the effects of the universal radiation field on the intensity 
of cosmic 7 -rays. The attenuation process of importance here is the pair- 
production process (53). This process can only take place if the total energy of 
the photons in the C.M.S. of the interaction is greater than or equal to 2mc 2 . 

The cross section for reaction (53) can be calculated using quantum electro- 
dynamics and a derivation may be found in Jauch and Rohrlich (1955). The im- 
portance of reaction (53) was first pointed out by Nikishov (1961) in considering 
interactions of 7-rays with ambient starlight photons and with the discovery of 
the universal radiation Gould and Schreder (1966, 1967a, b) and Jelly (1966) were 
quick to point out the capacity of the universe to y -rays of energy above 10 14 eV. 
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Stecker (1969b) and Fazio and Stecker (1970) generalized these calculations by 
including cosmological effects. 

In discussing reaction (53), we will generally follow the discussion of Gould 
and Schre'der (1967a, b) with one important difference. At the time Gould and 
Schre'der published their papers, it was generally thought that the universal 
radiation field had a somewhat higher temperature than the presently accepted 
2.7 K. Therefore, the Gould and Schre'der results have been corrected here to 
correspond to a 2.7 K radiation field. 

Denoting C.M.S. quantities by primes and noting that for reaction (53) 

E' + = E' = E' ( 54 ) 

+ e 

e e 

we can determine the threshold energy for the reaction by noting that the 
relativistic four -momentum invariance condition in this case reduces to the 
form 

(2E;) 2 = (E 2 7 + E y ) 2 - I (P 2 . 7 c + p y c)| 2 

= (E 2 -7 + 2E 2 7 E y + E 2 ) - (E 2 7 + E 2 - 2E 2-7 E r cos 9 ) (55) 

= 2E 2 7 E y (l - cos 9 ) 

At threshhold, both the electron and positron are produced at rest in the 
C.M.S. of the interaction. The minimum energy required corresponds to a 
head-on collision (cos 6 = -1). Equation (55) then reduces to the relation 
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( 56 ) 


If we consider a typical blackbody photon to have an energy of approximately 
10 9 MeV, then from (56) we find a threshhold energy of approximately 2.5 X 10 8 
MeV for reaction (53). However, this threshhold is somewhat blurred due to the 
fact that the blackbody photons are not all at the same energy but have a Bose- 
Einstein distribution given by the well-known relation 


o(E 2 7 ) = 


7T 2 h 3 C 3 J _ e E 2.7 /kT 


and also various possible values of cos 8 must be allowed for. 

The cross section for reaction (53) is given by (Jauch and Rohrlich, 1955) 


\ m c* 


(1-/3; 2 ) (3 -A; 4 ) In 2/3' e (2 - 


where the C.M.S. velocity of the electron (positron) is given by 


(i‘ = 1 

e V 


The absorption coefficient is then 


sin 6 


(b r ) = JJ dE 2> 7 d& — - — - n(E 2 7 )cr(E 2 7 . E y ) (1 - cos 8) (60) 


Gould and Schre/ier (1967a, b) have reduced equation (60) to the form 
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where 


v = ■ 


(m e c 2 ) 

E r kT 


( 62 ) 


They find that the function f(v) has a maximum value ~ 1 at v - 1 and that 
f(v) has the asymptotic forms given by 


f (V) -* v In y ' ^ ^ 1 for v « 1 


and 


(63) 


f(v) 


(?r - (* ♦£♦■■■) 


» i. 


They have also calculated x yy for y-ray interactions with various photon 
fields in interstellar space. The results of their numerical calculation are 
shown in Figure 7. 

B. Absorption by Interactions with Matter 

There are. two types of interactions of importance to consider here. The 
first is the Compton scattering interaction 


e + y - e + y 


(64) 


which we have discussed in connection with y-ray production. 

Electrons play the dominant role in the Compton scattering of y -rays. For 
y-rays of energy E y » me 2 , almost all of the energy of the y-ray is absorbed 
and then we can consider that the y-ray has disappeared. 
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Figure 7. Absorption probability per unit distance by y+ y- e + +e" as a function of 
photon energy for y-rays passing through various universal radiation fields. The 
contributions from the optical (0 X and 0 2 ) infrared (IR) 2.7 k blackbody microwave 
(M) and radio (R) fields are shown. Absorption at lower energies from X-rays is 
negligible. Oj represents the contribution from the light of population II stars; 0 
that from population I stars (from Gould and Schre'der 1967a, b); M is corrected to 
2.7 K (Stecker 1971a). 




In some cases it is useful to define an "absorption cross section", a , such 


that 




O' - 


(65) 


where AE y is the average amount of energy transferred from the gamma-ray to 
the electron. It is then found that (Heitler, 1954) with e = e y /m e c 2 


a =ecr for e«l 

a c 


( 66 ) 


and 


a ~a 
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In 

€ ~6 

In 

2e + I 
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for € » 1 


(67) 


The second type of gamma-ray absorption process in matter that we must 
consider involves the conversion of a gamma-ray into an electron-positron pair 
in the electrostatic field of a charged particle or nucleus. If we designate such 
a charge field by the symbol CF, such an interaction may be symbolically written 


as 


7 + CF -* e + + e" + CF 

The conversion interaction, or pair -production as it is usually called, has 
a cross section that involves an extra factor of the fine structure constant, 
a = e 2 /ftc, since it involves an intermediate interaction with an electrostatic 
field. At nonrelativistic energies, this cross section is a complicated function 
of energy which must be determined numerically (see Heitler (1954) for further 
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details) however a closed analytic approximation for this cross section may be 
given for energies greater than 1 MeV, which happily corresponds to the energy 
region where the pair -production cross section becomes more important than 
the Compton scattering cross section in determining the gamma-ray mass 
absorption coefficient for hydrogen gas. 

For pair -production in the field of a nucleus of atomic number Z, the cross 
section for reaction (68) is given by 


cr 

p 




In 2e - 



Z 2 


(69) 


for 1 « e « ol~ 1 Z~ 1/3 , 

which is the energy region where electron screening of the nuclear charge field 
may be neglected. The no-screening case, of course, also holds for an ionized 
gas (plasma). 

In the energy region where the complete-screening approximation is valid, 

Z 2 (70) 

for e » a _1 Z _1/3 . 

The threshold energy for pair production in the field of an atomic nucleus 
is, of course 2m e c 2 . In the case of pair production in the field of atomic electrons, 
the threshold energy for is 4m p c 2 . Above this energy, the pair production cross 
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section must be modified to include the additional contribution of the electrons 
and this may be done approximately by making the replacement 

Z 2 - Z 2 (l + < 71 > 

in equations (69) and (70), where the quantity £ varies from - 2.6 at E y - 6.5 MeV 
to - 1.2 at E y - 100 MeV. For gamma-ray energies above 200 MeV, £ - 1 and 
the pair production cross section has the approximately constant value of 1.8 x 10" 26 
cm 2 , according to the results of Trower (1966). The values of the various cross 
sections for y -ray absorption in matter, as discussed in this section, are shown 
in Figure 8. 

4. y-Rays Observed from the Galaxy 

In this section, we will now turn our attention to the interpretation of pres- 
ently existing observations of cosmic y -radiation. Because of lack of time, we 
will restrict ourselves to two important topics, galactic diffuse y-radiation and 
the diffuse y -ray background. This means we a omitting a discussion of such 
important topics as 1. the crab nebula and other possible point sources of y -rays, 

2. positron annihilation in the galaxy and nuclear de-excitation y ^ray lines, 

3. cosmic soft y -ray bursts, 4. y -rays from supernovae, and 5. solar y-rays. 

Some of these important topics will be touched upon here by Professor 
Pinkau, Professor Colgate and Professor Reeves. I would also like to recommend 
the discussion of some of these topics in the reviews published in the Proceedings 
of the International Symposium and Workshop on Gamma Ray Astrophysics held 
at NASA Goddard Space Flight Center last May (Stecker and Trombka 1973). 
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COMPTON AND PAIR PRODUCTION CROSS SECTIONS (lO* 6 CM* 2 / ELECTRON) 
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Figure 8. Compton scattering (cr s ), Compton absorption (cr a ), pair -production (cr ) and total 
(o'a + CT P ) cross sections as a function of y-ray energy for absorption of y- rays in hydrogen 
gas (based in part on work of Nelmes, 1953) from Stecker (1971a). 



The galaxy is a roughly disk-shaped collection of about 10 11 stars and a 
lesser mass of gas and dust having a radial dimension of about 15 kpc. The 
half -width of the gas disk, according to 21cm radio measurements of atomic 
hydrogen, is about 100 pc with a general tendency for a somewhat larger thick- 
ness in the outer regions. The sun is located approximately 2/3 of the way out 
from the center at a radial distance of about 10 kpc. 

Apart from these general considerations, it should be kept in mind at all 
times that the galaxy is a very complex object whose detailed general features 
such as spiral arms have not been charted to undisputed accuracy. Indeed, the 
required accurate observations to accomplish this are in some situations ex- 
ceedingly difficult if not outright impossible. 

Radio astronomers have for years been mapping accumulations and features 
involving relatively high densities of atomic hydrogen in the galaxy. They do this 
by studying the 21cm emission line of atomic hydrogen caused by the hyperfine 
splitting of the l 2 S ground state of the hydrogen atom as a result of the inter- 
action between the magnetic moments of the proton and the electron in the atom. 
The separation between these levels is very small because the magnetic moment 
of the proton is almost 2000 times smaller than that of the electron owing to its 
proportionally larger mass. The frequency of the emission corresponding to this 
energy-level separation is in the radio range at 1420 MHz. 

Radio source intensities are usually given in terms of brightness temperature 
T (Kelvin) corresponding to the emission temperature of a blackbody on the 
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Rayleigh -Jeans (low frequency) side of the Planck distribution. Thus 

T(y) = IM < P2) 

b 2^k 

Sometimes these intensities are given in terms of antenna temperature 
which is equal to the brightness temperature multiplied by antenna efficiency. 

For a gas radiating at constant temperature T but with opacity k , the brightness 
temperature measured along the line-of-sight in a given direction is 

T = J d S Te" Jo ' <ds, = j ds Te‘ T < s > = T(1 - e‘ T ) (73) 

b J 

Thus, at small optical depths r « 1, 

T b 2? Tr (74) 

and at large optical depths r >> 1, 

y y (75) 

“ 1 

The opacity « can be expressed as the atomic opacity a times the atomic density 
n of absorbing atoms. 

Consider an atom with two state levels, the upper level u and a lower level 
Then according to Boltzmann's law the proportion in each level as a function of 
T is 

n u _ -hv Q /kT (76) 

e 
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where v Q is the transition frequency and ^ and g Jt are the statistical weights 
of the two levels. Thus, for the 21 cm transition where g u = 3g^ (triplet to 
singlet) 


n u „ "h^/kT 


_S 3e 


The total absorption coefficient 


(77) 
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an l 



(78) 


where the second term takes account of a stimulated emission. From equation 
(76) this reduces to 


- an ^ (1 - e 


-h v^/kT 


hv Q hv 0 

) - arte for « 1 

J t kT kT 


(79) 


Thus, in the optically thin case, from equation (74). 

hv 0 

T. = Tr = a — — f n pds 

b k J l 

which is proportional to the total number of atoms along the line-of -sight. In 
the case of 21 cm emission 

T b (21 cm) =: 5.49 x 10 -19 / n^ds (80) 

In general, however, it should be kept in mind that equation (80), which is valid 
for an optically thin gas, in general only provides a lower limit to the amount of 
gas along the line-of-sight. 
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A 21 cm map of the atomic hydrogen density as constructed by Garmire and 
Kraushaar (1964) and is shown in Figure 9. 

In general, line emission is broadened in frequency by the Doppler shift 
effect of gas moving with various radial velocities v y - fi y c along the line of 
sight. Then 

v = v 0 (l + /y < 82 > 

What is usually measured is a velocity profile determined by the frequency 
distribution of the line emission f(v' ) so that a distribution in brightness tem- 
perature versus velocity or frequency is constructed. Then from equation (81) 

N h = J n^ds = 1.82 x 10 18 /T b (v) dv < 83 ) 

From the two dimensional distribution n (£, v) where 4 is galactic longitude, 
one can use a rotation model of the galaxy to construct a positional map of dense 
regions of atomic hydrogen. Such a map is shown in Figure 10 with some of the 
21 cm atomic hydrogen "arms" of the galaxy named. 

One can also plot a graph of the average hydrogen density seen in 21 cm 
emission as a function of galacto -centric distance in kpc. Such a plot is shown 
in Figure 11. 

The galactic distribution of y-ray emission was first mapped by Kraushaar, 
et al. (1972) on OSO-3 with rather limited resolution. More recently, the galactic 
distribution was remeasured with better angular resolution by Kniffen, et al. 
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Figure 9. Contour map of N h> the number of atomic hydrogen atoms per cm 2 along the 
line of sight in galactocentric coordinates (from Garmire and Kraushaar (1964)). 




Figure 10. A smoothed spatial diagram of the locations of the maxima of the 
matter density deduced from the 21 -cm neutral H line measurements and the 
density-wave theory by Simonson (preprint). 
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(1973) on SAS-2. The results of both measurements are shown in Figure 12. 

The map of Kniffen, et al. has recently been updated by the addition of new data 
having been analyzed and the updated distribution is shown in Figure 13 (Fichtel 
1974). 

The data from SAS-2 of y-rays above 100 MeV shows radiation coming from 
the galactic disk. The galactic latitude distribution of the emission is shown in 
Figure 14. Figure 14 indicates that the radiation in the anticenter direction 
comes from within 6° of the galactic plane whereas that in the direction of the 
galactic center also is limited to within 6° with perhaps a particularly intense 
component restricted to within 3° of the plane. Since 3° was the angular resolu- 
tion of the instrument, this narrow source could have a true width of less than 3° 
corresponding to a source distance of at least 2 to 4 kpc from the sun toward the 
inner galaxy. A two dimensional SAS-2 map of the galactic disk is shown in 
Figure 15. 

Data obtained by Samimi, et al. (1974) on the latitude distribution of the 
galactic y -radiation also indicate a width of 3° for most of the radiation. 

The longitude data in Figure 13 show a broad flat region of intense emission 
within 30° to 40° of the galactic center on either side. This indicates that there 
is a large emission rate within 5 to 6 kpc of the galactic center. 

Puget and Stecker (1974) have geometrically unfolded the SAS-2 distribution 
shown in Figure 12 and the resultant emissivity as a function of galactocentric 
distance is shown in Figure 16. Strong (1974) has performed a corresponding 
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GALACTIC LINE FLUX (ARBITRARY UNITS) 


C» 



GALACTIC LONGITUDE 


Figure 12. Distribution in galactic longitude of the galactic y-ray line flux as 
observed by OSO-3 (Kraushaar, et al. 1972) and SAS-2 (Kniffen, et al. 1973) 
normalized arbitrarily for purposes of comparison. Statistical uncertainties 
quoted for these results (not plotted) are significantly greater for the OSO-3 
results as compared with the SAS-2 results. The interpolation of the SAS-2 
data was obtained on the assumption of a smooth average variation in the 
y-ray production rate (Puget and Stecker 1974). 
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RELATIVE NUMBER 



Figure 13. Distribution of high-energy (>100 MeV) gamma -rays along the galactic 
plane. The diffuse background level is shown by a dashed line. The SAS-2 data 
are summed over |b n | < 10°. The ordinate scale is approximately in units of 10 
photons cm -2 rad -1 sec -1 (Thompson, et al. 1974). 





Figure 14a. Distribution of high-energy (> 100 MeV) gamma -rays summed over 335° < l n < 25° 
as a function of b 11 . The solid line represents the sum of the two distributions with equal areas, 
one representing only the detector resolution, the other a gaussian with a width o- = 6° {Thompson, 
et al. 1974). (b) Distribution of > 100 Mev gamma-rays summed from 90°<^ II < 170° and 200° 
< ^11 < 270° where data exists, (Thompson, et al. 1974). \ 



DECLINATION 



Figure 15. Schematic diagram of relative intensity of the gamma radiation 
above 100 MeV as deduced from preliminary data from SAS-2. Increasing 
shading indicates higher intensity, but no single (4°x4°) box has sufficient 
numbers of gamma rays to justify statistically conclusions to be deduced 
for it alone. It is clear, however, that the galactic plane stands out sharply 
(Fichtel, Hartman, Kniffen, and Thompson, 1973*). 

"Gamma Ray Astronomy in the Time of SAS-2”, invited talk 140 AAS Meeting 
Columbus, Ohio, June 1973. 
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Figure 16. The value for hQ (S) (7-ray emissivity times disk width) 
given by Puget and Stecker (1974). A constant value for hQ is as- 
sumed for R G Si > 8 kpc. Uncertainties in the determination of the 
value of hQ(S) as evaluated from an integral equation of Puget and 
Stecker grow quite large in the central region R 0 Si < 4 kpc. 
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unfolding of combined data from SAS-2 and OSO-3. He obtained a similar distribu- 
tion to that shown in Figure 16. The results show a maximum in the y-ray 
emissivity at ~5 kpc and a reduced emissivity within 3 kpc of the galactic center. 
An unfolding of Figure 13 would probably show a less pronounced maximum at 
5 kpc, but this unfolding has not been attempted here because the SAS-2 group 
will be publishing their final and complete longitude data in the near future. 

Data have also been obtained on the energy spectrum of the galactic 
y -radiation in the direction of maximum intensity by Kniffen, et al. (1973). 

These spectral measurements on the integral y-ray spectrum are shown in 
Figure 17 to be consistent with a two-component origin with -^70 percent of the 
radiation above 100 MeV due to 7r°-decay as shown in figure 4 and ~ 30 percent 
due to Compton radiation from cosmic ray electrons interacting with the higher 
radiation field in the inner galaxy (Stecker, et al. 1974). Under this interpretation, 
one would expect that almost all of they -radiation above 100 MeV in the. outer 
regions of the galaxy would be due to 7r°-decay. This interpretation would also 
be consistent with the y-ray production rate deduced by Kraushaar, et al. of 
(1.6 ± 0.5) x 10 -2S cm' 3 s” 1 above 100 MeV which is only slightly larger than the 
theoretical calculated value of (1.3 ± 0.2) x 10 -2S n H cm -3 s" 1 for the solar vicinity 
(Stecker 1973) with ~ 0.7cm~ 3 for atomic hydrogen alone. 

While this interpretation presents a coherent and plausible model for dis- 
cussion, one should bear in mind that the empirical situation is still in a state 
of flux and one measure of uncertainty is indicated by the various different 
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Figure 17. Comparison of SAS-2 (Kniffen, et al. 1973) spectral data on 
7-radiation from the inner galaxy with a two-component model based 
on 70 percent pion decay (Stecker 1970) and E" 1 integral Compton spec- 
trum (Stecker, et al. 1974). 
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measurements of the inner-galactic flux as shown in Figure 18 and the differ- 
ential data of Samimi, et al. (1974) as shown in Figure 19. In particular, the 
data summarized in Figures 18 and 19 appear to support the presence of a soft 
component of radiation at low energies possibly due to Compton interactions or 
bremsstrahlung. 

It is instructive to estimate Compton and bremsstrahlung intensities in the 
galaxy relative to 77°-decay y -radiation. First a crude estimate of the Compton 
radiation in the galactic disk. Starlight photons in the galaxy have an average 
energy of ~ leV and a number density of ^0.4cm'3 (Allen 1973). Let us consider 
the total gas density of atomic and molecular hydrogen to be ~1 cm" 3 . As pointed 
out previously, the mean energy of cosmic ray protons which produce 7r° -mesons 
leading to y-rays in the observed energy range £ lOOMeV is ~3GeV. By a 
curious coincidence, a 3 GeV electron with y ~ 6 x 10 3 interacting with a leV 
photon will produce a y -ray of energy 

E y =i y 2 eV - 50 MeV (84) . 

o 

so that electrons of comparable (or slightly higher) energies will produce y-rays 
in the same energy range as protons which produce tt ° decay y -rays. The ratio 
of electrons to protons at these energies is ~ 10‘ 2 . The ratio of total y- ray inten- 
sities is then 
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Figure 18. Summary of the integral flux measurements for the galactic 
center region given by various experimental groups. 
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Figure 19. The differential photon spectrum of y-rays as 
determined from the emulsion measurements of Samimi, 
et al. (1974) with the detector response unfolded, a) only 
y rays coming from the enhanced 3° band -3° < b :i < 0°, 
b) atmospheric y-rays. 
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(85) 


More involved calculations using equations (11) and (12) give similar results. 
Equation (11) and (12) can also be used to calculate the y-ray spectrum produced 
by interactions of electrons with the 2.7K universal background radiation. The 
ratio of 2.7K blackbody Compton y -rays to starlight Compton y -rays is 


T /e- \r-3/2 . . 

*2. 7, c _ ^2-7 | E 2 ■ 7 \ ^ f o. 2 5\ 

I., e " P* ) " V 0.4 7 


(6.4 x 10 ‘ 4 ) -0 ' 2 2? 2.7 


( 86 ) 


As can be seen from Figure 13, they -ray intensity in the inner galaxy is 
about 5 times that in the anticenter direction when one subtracts out the peaks 
due to the crab nebula at 180° and the Vela source (pulsar?) at 270°. A similar 
large ratio was found by Kraushaar, et al. (1972) on OSO-3. The OSO-3 ratio 
was speculated by Steeher and Stecker (1970) to be due to two effects: 1) 77 0 
production off H 2 not seen in 21 cm emission and 2) a Compton source from a 
large infrared radiation field at the galactic center (Hoffman and Frederick 1969). 
However, the Compton source would produce a clear peak at -l - 0° whereas the 
new SAS-2 data do not indicate such a peak but rather a broad flat region of 
relatively intense radiation (see Figure 13). 
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Bignami and Fichtel (1974) have attempted to explain the SAS-2 longitude 
distribution as due to an n 2 enhancement of the product of gas density n and 
cosmic -ray intensity I with large density contrasts in galactic spiral arms as 
defined by 21 cm measurements of atomic hydrogen. Their results are shown 
in Figure 20. 

Stecker, et al. (1974) and Puget and Stecker (1974) pointed out that the ap- 
parent peak in y-ray emissivity at ^5kpc implied by the data of Kniffen, et al. 
corresponded to a maximal dissipation of the kinetic energy of outward gas motior 
in the galaxy and that possible cosmic-ray compression in that region and first- 
order Fermi acceleration could enhance the cosmic-ray density in that region to 
a large enough extent to explain the enhancement in y-ray emissivity. 

An enhancement of cosmic-rays in a localized region of the galaxy can be 
due to three factors (1) an increase in the density of cosmic-ray sources (or 
alternatively the production rate) in the region (2) an increase in the trapping 
time (escape time) of cosmic-rays in the region, and (3) acceleration and com- 
pression of cosmic-rays in the region. We may write this as follows: 

ln (I CR /I 0 ) - ln(Q CR /Q a ) + ln(T CR /T 0 ) + S (87) 

where the first term on the right hand side of equation (87) represents the en- 
hancement in the production rate, the second term represents the trapping factor 
and the third term represents the effect of acceleration and compression. The 
first term can only be guessed at, given our present lack of knowledge of the 


56 



GAMMA RAYS (E>IOO MeV)/(CM -SEC- SR) 



Figure 20. Longitudinal distributions of galactic gamma-flux integrated over ±10° in 
b 11 from Bignami and Fichte! (1974). SAS-II points are given together with their 
error bars (Kniffen, et al. 1973). The thick line represents the model of Bignami 
and Fichtel smoothed in 10° -fc 11 intervals. The thin line represents the model in 2° 

intervals. The dash line ( -) gives the contribution of the Sagittarius and Norma- 

Scutum arms and dash-dot (-.-.-.) the contribution of the Sagittarius arm alone. 




ultimate origin of cosmic rays. However, as an indication of the possible effect 
of source density, one can note the distribution of supernova remnants in the 
galaxy which may be proportional to the density of cosmic -ray sources if we 
assume that cosmic-rays are produced by supernova explosions (Ginzburg and 
Syrovatsky 1964) or in remnant pulsars (Gunn and Ostriker 1969). According to 
Ilovaisky and Lequeux (1972) the density of supernova remnants increases roughly 
by a factor of 2 over the local galactic value for w < 0.8 and drops sharply for 
£ > 1.2. A more recent study by Clark, et al. (1973) is probably less susceptible 
to selection effects because their survey included remnants down to lower lum- 
inosity levels. The results obtained by Clark, et al. confirm the earlier results 
of Ilovaisky and Lequeux regarding the general distribution of supernova remnants 
in the galaxy. We will, therefore, estimate here thatQ cr /Q 0 ~ 2 for the enhanced 
region 0.4 < a> <0.6. This factor, by itself, cannot account for the order-of- 
magnitude enhancement deduced for I /I . 

The second factor, T cr /T 0 is so difficult to estimate that we will treet it as 
a free parameter > 1 to be solved for. It is not unreasonable to expect more 
effective trapping in the inner galaxy (i.e. T cr /T 0 >1) due to compression and a 
resultant stronger magnetic field strength. 

The factor S can be broken up into two parts, acceleration factor S ACC 
which accounts for acceleration of the more numerous lower energy cosmic -rays 
(given an assumed power-law differential energy spectrum of the form * E~J" ) 
to an energy above the threshold for production, and a simple density 
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enhancement S D due to the lower volume of the compressed retion in which the 
trapped particles find themselves. If we designate the specific volume compres- 
sion rate by 

as- (1/V) (dV/dt) < 88 > 

then obviously 

§ = aT (89) 

D CT 

The acceleration factor can be estimated therm adynamic ally. Regardless 
of the exact details of this process, which may be coherent, first order Fermi 
acceleration which can transfer the momentum of moving "clouds" with trapped 
magnetic irregularities to cosmic— rays trapped in the compressed region (Fermi 
1954, Stecker et al. 1974). The increase in the energy of the individual cosmic- 
rays can be treated as an adiabatic compression heating of a "cosmic-ray gas 
molecule". The energy enhancement factor as 

ln(E/E 0 ) = (7 - 1) aT cr 

where y = 4/3 for relativistic cosmic-rays y = 5/3 for sub-relativistic 
rays. For a cosmic -ray differential energy spectrum of the form I £r a 
then find 

S ACC = < r - 1 )^- 1 ) aT cr 

and therefore (®1) 

S = [<T - 1) (7 - 1) + 1] o.T ct 


(90) 

cosmic- 

-P 

E , we 
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(Stecker, et al. 1974). The value of the specific volume compression rate a(oJ) 
can be given in terms of the radial expansion velocity V f of the gas deduced 
from the 21 cm observations of the "expanding arm" feature. Under this 
interpretation 


R 0 " = 


dv 

r 

d£ 


(92) 


where the galactocentric distance R = 

The function a (25) deduced from the observations of v (23) (Shane 1972, Sanders 
and Wrixon 1972, 1973) is shown in Figure 21. It is positive and maximal in the 
region of observed maximum 7 -ray emission (see Figure 16) and is negative in 
the inner region 23 < 3 kpc where there may be a significant drop in 7 -ray emis- 
sion, although the uncertainties in that region are very large. 

Because of the large uncertainties in T cr and T 0 we will consider a range 
of values for T cr and T cr /T Q . We will consider 3 x 10 6 < T < 3 x 10 7 yr 
where the upper limit on T cr is taken to be the deduced age of the expanding 
feature (Oort 1970, Van der Kruit 1971). The cosmic-ray lifetime in the local 
region, T s , is given by O'Dell, et al. (1973) and Brown et al. (1973), to be 10 6 yr < 
T cr £l ° 7 yr. Takinga= 2 x 10 ~ 15 s " 1 in the region of maximum compression (see 
Figure 8 ) and assuming I cr /I 0 hr 15 and Q cr /Q a - 2 , P = 2.5, S = 1.75 a T cr = 
1.5 x lO" 15 T cr , the value given for the parameters T /T_, e s and T and T 

cr u cr © 

consistent with equation (87) are given in Table 1. The factor e 5 represents the 
enhancement effect due to acceleration and compression. 
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Table 1 


Lifetime of Cosmic -Rays in the 4 to 5 kpc Region Calculated 
from the Acceleration-Compression Model 


T /T n 

e s 

T cr (10 6 y) 

% (io 6 y) 

2 

3.7 

12.6 

6.3 

3 

2.5 

8.7 

2.9 

4 

1.9 

6.0 

1.5 


Aside from the possibility that cosmic-ray enhancement accounts for the 
large 7 -ray emissivity in the 5 kpc region, recent carbon monoxide emission 
measurements in a 2.6 mm microwave survey of the galaxy by Scoville and 
Solomon have given new life and strong emphasis to the idea that much if not 
all of the increase in the 7 -ray emmisivity could be due to interaction between 
cosmic rays and large amounts of H 2 not seen in 21 cm emission (Stecker 1969c, 
1971; Stecher and Stecker 1970). However, the recent measurements of Scoville 
and Solomon point up this possibility in a much more specific and meaningful 
way, since their measurements indicate an increased enhancement in CO and by 
implication an increased density in the very region where the y-ray emis- 
sivity appears to be enhanced. This has led Solomon and Stecker (1974) to 
suggest that both enhancements provide evidence for the existence of a galactic 
H 2 "ring" or "arm" feature at - 5kpc. 

Molecular hydrogen is expected to be the predominant form of hydrogen in 
cool clouds of sufficient density (Solomon and Wickramasinghe (1969), Hollenbach 
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and Salpeter (1971), Hollenbach, Werner and Salpeter (1971). However, it is 
difficult to measure its galactic distribution directly. Strong H 2 absorption lines 
have been seen in the UV in almost all nearby clouds by the Copernicus satellite 
(Spitzer, et al 1973) but UV observations of H 2 at distances greater than 1 kpc 
from us are not feasible because of the large extinction of UV radiation by inter- 
stellar dust. 

Because H 2 has no permanent dipole moment, it cannot emit electric dipole 
radiation and, even though quadropole vibration-rotation features can in principle 
be observed in the infrared, such features are inherently very weak. 

CO emission can be used as a tracer of H 2 in molecular clouds because the 
most important source of CO excitation in these clouds is by collisions with H 2 . 

The galactic longitude distribution of CO emission in the galactic plane 
obtained by Scoville and Solomon (to be published) is shown in Figure 22. The 
data represent the antenna temperature integrated over the velocity profile of 
the 2.54mm CO emission line and would therefore be proportional to the number 
of CO molecules along the line-of-sight in the optically thin case as we discussed 
previously with the 21 cm line. Measurements were also obtained of the velocity 
profile function T a (-£, v) for each 1° of longitude from 0° to 90° and sample profiles 
are shown in Figure 23. 

The function T A (£ , v) can be converted into a galactocentric distance distri- 
bution by using a model for galactic rotation (Schmidt 1965). Such an unfolding 
is shown in Figure 24 which can be compared with the y-ray emissivity function 
given in Figure 16. 
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Figure 22. The intensity distribution of 2.6mm line emission in the galactic plane from the 
J = 1 - 0 transition of carbon monoxide integrated over velocity as a function of galactic 
longitude (Scoville and Solomon, to be published; Scoville, Solomon and Jefferts 1974). 
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Figure 23. Sample CO line emission spectra obtained in the survey of 
Scoville and Solomon (to be published). 
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Solomon and Stecker (1974) have used the CO emission data and the y -ray 
emission data to place limits on the mean density of H 2 in the 5 kpc feature. 
They have obtained 

1 < (5 kpc) < 5 molecules/cm 3 (^3) 

The lower limit closely corresponds to the optically thin case for CO emis- 
sion and the upper limit corresponds to optically thick emission of the type ob- 
served in nearby dark clouds. In the optically thick case for C 12 O in nearby 
clouds, emission from C 13 0 is optically thin and one can use the abundance ratio 
of C 13 /C 12 to deduce that C 12 O is in the range 10-20. It should be kept in mind 
that in order to obtain enough collisional excitation of CO to the J = 1 rotational 
state to produce a measurable amount of J = 1 - 0 2.64 mm emission, the H 2 
density in the cloud must be in the range 10 3 < H 2 < 10 4 cm 3 . Therefore, 
equation ( 93 ) represents only a smoothed out mean of a very uneven distribution 
of H 2 . 

It can thus be seen how 7 -ray astronomy gives us another window on the 
study of galactic structure and dynamics. The various implications of the SAS-2 
data we have just discussed represent only a first step toward our understanding 
of 7 -ray emission in the galaxy. Improved y-ray telescopes will someday pro- 
vide us with a clearer picture of the galaxy. 
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5 . Extr agalactic y -Rajs 

A. Redshifts and Cosmology 

Much data have recently been obtained on the diffuse y - ray background radi- 
ation which appears to be isotropic in origin. These data have now defined a 
continuous background spectrum up to an energy of 200 meV. They are sum- 
marized in Figure 25.* 

The cosmological nature of this radiation leads us to examine its origin in 
the context of the expanding universe model where such radiation is redshifted to 
lower energies as we see it because of the Doppler effect. 

The simplest expanding universe model is one in which the universe is both 
homogeneous and isotropic. Such a model is quite adequate for our purposes here 
and support for it comes from the remarkable isotropy of the 2.7K microwave 
background radiation which is believed to have originated much earlier in the 
history of the universe than the observed y -ray background. A homogeneous- 
isotropic universe can be described by the Robertson-Walker space-time metric 
of the form 

ds 2 = c 2 dt 2 - R 2 (t) du 2 (94) 

where R(t) is a scale factor which describes the expansion of the universe as a 
function of time according to the solution of the Einstein equations of general 
relativity. 

* 

Recent data reported by Tanaka (1974) in the energy range up to 7.5 MeV are not shown in the 
figure but are consistent with other data in this range. 
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Figure 25. Observational data on the gamma-ray 
background energy spectrum. The highest energy 
point of Vette, et al. (1970), shown with a dashed 
line, and possibly the neighboring point are now 
thought to be erroneously high due to an ineffi- 
ciency in the anticoincidence circuit of their de- 
tector which should not significantly affect the 
points at lower energies (Vette, private communi- 
cation). 
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Photons travel along the null geodetic which obeys the relation ds 2 = 0, i.e., 


from (94) 



(95) 


with t e being the time the photon was emitted and t r being the time the photon is 
received. 

The emitting and receiving points are embedded in the metric so that the 
distance between them is changing by the scale factor R(t); the dimensionless 
metric distance, u, is a constant. If we therefore consider two successive wave 
crests of a light ray as being emitted at times t e and t e +A t respectively and 
being received at times t f and t f + At , then 


Thus 



t +At 
r r dt 


t +At 

e e 


R(t) 


u = const. 


(96) 



+At 

At 


e 


dt 

R(t) 





dt 

R(t) 


(97) 

At At 

— - A = 0 

R(t r ) R(t e ) 


At At 

e _ r 

R(t7) = RCT) 
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Since the wavelength of the emitted wave is cAt e and that of the wave when 
received is cAt r , the wavelength is shifted by the amount 


or 


_A\ R(t r ) - R(t e ) 

\ R(t e > 


( 98 ) 


R(t r ) 

R(t e ) 


1 + z 


We have observed this shift in the spectral lines of distant galaxies as always 
being toward longer wavelength so that R(t r ) > R(t e ). From this evidence, it 
has therefore been deduced that our universe is expanding with time. 

B. Gamma- Ray Fluxes 

Let us now consider the effect of cosmological factors in calculating gamma- 
ray fluxes emitted at large redshifts, z. The number of photons received per 
second is reduced by a factor R(t e )/R(t r ) from the number produced per second 
at time, t e . We consider here gamma rays produced in particle collisions between 
two components having densities n a (t e ) and n b (t e ) respectively. We specify the 
differential photon intensity produced per collision as 


G(E y ) (cm 2 • s • Sr -MeV- cm -6 ) -1 

Then the differential photon flux received at t r is given by 

47m (t e ) n^t.) G(E y e ) dE y e dV e dt e 

dF = — 

47rR 2 (t r )u 2 


( 99 ) 
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where the numerator represents the photon flux emitted at t e , and the denominator 
indicates the fact that at t r this flux is evenly distributed over a spherical wave- 
front of radius R(t r ). We now define the three dimensional length, 

d l = R(t) du (100) 

so that the volume element 

dV e = d-£-[R 2 (t e ) u 2 dO] ( 101 > 

Since 

dt e = CR(t e )/R(t r )] dt r 

and (102) 

dE r >e = [R(t r )/R(t e )] dE y r 

because the energy of a gamma-ray is inversely proportional to its wavelength, 
we may substitute (101) and (102) into (100) and obtain 

dF r = n a(t e ) "bCU G{[R(t r )/R(t e )] E r>r } 

(103) 

47rR 2 (t ) uMfid'tdE dt 

x v e y y, r r 

47rR 2 (t r ) u 2 

By making use of equation (98) and dropping the subscript, r, since we only 
measure gamma-rays when they are received, equation (103) reduces to 
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dF 


dfidtdE 


dl = 


7 


n a (z) I\(z) G [(1 + z) dl 
(1 + z) 2 


(104) 


n a (z) n„(z) G[(l + z) E y ] 


(1 + z ) 3 


( 

\dz/ 


dz 


Equation (104) is quite useful in evaluating the metagalactic gamma-ray spectra 
from various high-energy interactions. The results are obtained from numerical 
integration of the relation 


I(E„) 


-f 


dz n (z) n,(z) 


G[(l + z)E y ] Atf 


/ d^\ 

\dz/ 


(105) 


(1 + z ) 2 

where the factor dl/dz is determined from the Einstein field equations to be 


dl C ,, .-2 n 

— =— (1 + z) 2 (1 + Oz) 

dz H 0 


(106) 


where H 0 is the Hubble constant and 


Q= — 

n 


(107) 


where n 0 is the present mean density of all the matter in the universe and n c 
is the critical density needed to gravitationally close the universe. 


The value for cH' 1 is in the range 


10 28 < cH; 1 < 2 X 10 28 cm 


(108) 


and 


3 x 10~ 6 < n < 10 -5 cm* 3 


*\j C 


(109) 
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If absorption effects are included, equation (105) becomes 


I(E r ) 


c r - * , x g[ ( i+z )v -t(e o 

= dzn a (z)n b (z) e 

n o •'o 


(HO) 


(1 + z) 4 /l + Qz 

Since the universal blaekbody radiation plays an important role in various 
high-energy production processes, it is important to determine how this radiation 
behaves as a function of redshift. At the earliest stages in the evolution of the 
universe, when the matter in the universe was in thermal equilibrium with the 
universal blaekbody radiation owing to Compton interactions between thermal 
photons and electrons, T m was equal to T r (the temperature of the radiation 
field). Zeldovich, et al. (1969) have shown that even though the intergalactic gas 
has cooled to the 50% neutral point by the time corresponding to a redshift 
z = 1200-1300, the small fraction of ionized material left at lower redshifts is 
enough to sustain temperature equilibrium between matter and radiation field 
until a redshift of 150-200. 

The photon density of a radiation field at temperature T 0 is given from the 
Planck formula as 


n r.o( £ ’ V de = 


877 

h 3 c 3 


e 2 de 


e/kT 


( 111 ) 


r, 0 


- 1 


At a redshift z \ 0, this distribution is then given by 


n (z, £) de = JH J 1 + z > 3 e2dg (112) 

h 3 c 3 e (i+«)e/kT y>0 _ 1 

= (i + = 0 3 n r , 0 [ 6 > T r ( z ) de 
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where 


T r (z) - (1 + z) T r o 


(113) 


so we find that the redshifted Planckian maintains its form with the parameter 
T(z) in the exponent given by equation (113). The photon number density redshifts 
in the same manner as the matter density viz. n “ R(t) 1 « (1 + z) . The 
energy density in the radiation field then redshifts as 

P t = P r,<>( 1 + Z > 4 < 114 > 

in accord with the definition of T r (z) given by equation (11-12) and the relation 
p * T 4 . 

We therefore find that 

T ( Z ) = T (z) = T r p (l + z) (115) 


for z > 150 - 200 


At lower redshifts, when the matter has thermally decoupled from the radi- 
ation field, the momentum distribution of the atoms in the intergalactic gas is 
given by the Maxwellian, 


n Q (p, T m,o) d P = COnSt - >< P 2e 


-p i /2mkT_ 


dp 


(116) 


These atoms lose momentum through collisions due to the effect of the 


overall adiabatic expansion of the universe; their resultant momentum change is 
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P(Z) = Po R(I) = Po(1 + Z) 


(H7) 


Substituting (117) into (116), we then find 

n ro (z> P) dp = (1 + z) 3 n m 0 [p, T m (z )] dp (118) 

with 

W = VoC 1 ^) 2 ( 119 ) 

From equation (118), it is immediately evident that the total number of 
atoms, ~nR 3 is conserved during the expansion as must be the case. 

We thus find that 

T m (z) = T ro o(i + Z ) 2 for z< 150 - 200 (120) 

with 

T m (150 to 200) - T(150 to 200) ( 121 ) 

The adiabatic expansion concept may be used to link the derivations of 
equations (113) and (119) for the temperature -redshift relations. In an adiabatic 
expansion, it is well known that 

TV^ _1 > = const. (122) 

where the quantity, y , here is the ratio of specific heat at constant pressure to 
that at constant temperature and V is the volume of the gas (in this case the 
volume of the universe) which from (98) is given by 

V(z) = V 0 (l + z)" 3 ( 123 ) 
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For photons andmonatomic gases at relativistic temperatures (pc - E), 
y - 4/3, whereas for nonrelativistic gases y = 5/3. Thus, from (123) and (122) 
one can immediately obtain both equation (113) for photons and relativistic gases 
and (119) for nonrelativistic gases. 

We can use these equations to take into account the cosmological effects in 
computing the extragalactic gamma-ray spectrum from Compton interactions 
k e tw ee n cosmic-ray electrons and photons of the universal blackbody radiation 
field in intergalactic space. The local (z = 0) spectrum is given from equation 
(11) as 


Vo (e y ) =| ff, cH ; lf (T) (me 2 )' 1 ' 1 "’ (3.6k) 


fr) 


^e.O^r.O 


(124) 




, o y 


for Compton gamma-rays produced by cosmic-ray electrons having a power law 


energy spectrum of the form 


I(E e ) = K ei0 E 


-r 


(125) 


Under the conditions where (125) is valid over all redshifts from 0 to z 
and where electron energy losses other than those from the universal expansion 
are neglected, the transformation from K e 0 to K e (z) is given by 
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(126) 


K e<» = Ke.oC 1 + Z > 3 C 1 + zf ' 1 


= Ke.oC 1 + z ) 


r+a 


where the first factor of (1 + z) 3 represents the density effect and the factor of 
(1 + z f~ l represents the transformation of the power-law energy spectrum from 
a burst of relativistic cosmic -ray electrons occurring at z max > z. 

Using equations (114) and (115), as well as (126) and substituting them into 
the general formula (105), we find 

h&y) = K.o&y) f maX dz (127) 

Jq V'l + Qz 

so that under these conditions, the power-law form for the Compton y-ray 
spectrum is maintained but the cosmological flux is enhanced by a factor given 
by the integral over z in equal (124). 

However, Brecher and Morrison (1967) showed that the true situation is not 
as simple as that given by equation (127) under the conditions when the high 
energy end of the electron spectrum is steepened owing to energy losses from 
the Compton interactions themselves. Since the blackbody radiation density is 
proportional to (1 + zf , electrons lose energy primarily from Compton inter- 
actions with the energy loss rate being given by 



4 CT T C ^r 


y 


2 _ (1 + z ) 4 ^2 


m cr 


(128) 
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with 


T,, ■ 3m - c; = 7.7 x 1 0 19 s <129) 

4 ^T C P r ,0 

The energy loss from the universal expansion is given by 

(< ty \ _ y (130) 

\dz/ 1 + z 

exp 

From (130) and (106), we obtain 

/dy\ = /&A fdt\ = (1 + z) 2 y 2 (131) 

\dz) c \dt j c \dz J Hq t 0 v/1 + Qz 

The two energy loss rate terms are equal at the critical energy, E c = 7 C m e c 2 , 
where 

V; vl 1 E 256(1 t_n Z )' /2 (132) 

>C " (1 t z)3 (1 t z? 

or 

E ~ 130^1-ti^l MeV ( 133 ) 

(1 + z ) 3 

For E < E , the electron spectrum maintains its power-law form since 

C 

these electrons lose their energy through the universal expansion and equation 
(126) is therefore valid. However, for E > E c , the equilibrium electron flux is 
given by the solution of a rate balance equation under the conditions when Compton 
losses dominate. Under these conditions, 
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By 


(134) 



(1 + z) 4 y 2 ~ | 



n 

(Brecher and Morrison, 1967), where k q y is the original (injection) electron 
spectrum and P is the exponent of the resulting equilibrium electron spectrum. 
It follows from equation (134) that 


r = r + 1 ( 135 ) 

and 

K = kqT ° (136) 

P(1 + z) 4 

so that the equilibrium electron spectrum is depleted by a factor proportional to 
(1 + z) 4 by Compton interactions with the universal radiation field and, in addi- 
tion, the exponent of the electron spectrum is steepened by one power of E . 

This steepening in the electron spectrum corresponds to a change of 1/2 in 
the exponent of the Compton gamma-ray spectrum at an energy 


r.c 


— £, 


y 


2 _ 


33 


1 + fiz 
(1 + z ) 5 


eV 


(137) 


Thus, the Compton gamma-rays are expected to be produced by the steepened 
electron spectrum at all gamma-ray energies. From (136), Brecher and Morrison 


conclude that unless there is a large evolutionary factor in the electron production 
spectrum, i.e., k q ~ (1 + z) m , where m > 4, there will be no significant enhance- 
ment of Compton gamma-rays at large redshifts over those produced at the present 
epoch. 
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The local gamma-ray spectrum from relativistic bremsstrahlung inter- 
actions is given by (36) as 


I b,o( E r ) = 3 - 4x 10 " 26 n o cH o 


,-i M > E r > 


(138) 


Making use of equations (105) and (126), we obtain for the cosmological 


bremmstrahlung production spectrum 


^b^yO' ax) 


3.4 x 10“ 26 n Q cH^ 1 K 

r - i 


e , 0 


X 



(1 + z) 3 (1 + zf* 2 [(1 + z) E y0 l 
(1 + z ) 4 A 4 fiz 


3.4 x 10“ 3 Q 

m 


K 


e , 0 




dz 


(1 + 2 ) 

vTTHz 


(139) 


6.8 x 10“ 3 v ,,-r 
~Z ~ i N.O^y.O 


f i + £ _ i_] yi + n z _ (l - 

i 3 3fi/ [ 3C1J 


In particular, for the Einstein-de Sitter universe where fi = 1, 


W = I Ib . o ( Eyo ) [(1 + :,. x ) 1/2 - U 

and for the low density model where fi z « 1, 

IbCEyO' ^B*) = 2 ^b.o^yo) t^ 1 ' + Z max^ ~ ^ 


(140) 


(141) 
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Where electron steepening by Compton interactions with the 2.7K radiation 
is important, the situation is more complicated and a more detailed treatment of 
this problem has been given by Stecker and Morgan (1972). 

6. Gamma-Ray Absorption Processes at High Redshifts 

In Section 4, we discussed in detail the processes which result in the ab- 
sorption of cosmic gamma-rays. In that section, we pointed out that there are 
three main absorption processes of importance under the conditions in inter- 
stellar and intergalactic space: 1) pair -production interactions of cosmic 
gamma-rays with the universal blackbody radiation field, 2) Compton interactions 
of cosmic gamma-rays with electrons the interstellar or intergalactic gas, and 
3) pair -production interactions of cosmic gamma-rays with atoms of the inter- 
stellar or intergalactic gas. In this section, we will show how the effectiveness 
of these absorption processes is enhanced at high redshifts due to increased 
photon and matter densities at these redshifts when the universe was in a more 
compact state. We will also derive the energy-dependence of these absorption 
processes, taking cosmological factors into account. 

We begin our discussion by considering the absorption of cosmic gamma- 
rays by pair-production interactions with photons of the universal blackbody 
radiation field, viz., reaction (53) whose absorption coefficients are given by 
equations (61)-(63). 

For cosmological applications, we must take into account the redshift 
dependences of T and E y in an expanding universe, 
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T - T 0 (l + z) 


and (142) ' 

E y = E y0 (l + z) 

where the subscript zero refers to presently observed (z = 0) quantities, so that 
T 0 = 2.7 K. 

Taking the z -dependence into account, we then find that the condition v » \ 
in equation (63) is applicable in the energy range 


E « 
r 


1.12 x 10 6 GeV 


(143) 


(1 + z) 2 

The optical depth of the universe to gamma-rays is then given by 


r(E.., z ) 

v y max 7 


f 


max max ' 


AlK yy & r z) 


(144) 


■f 


dzK yy (Z y , 


z) 



where, from equation (106) 


d'C ^ 10 28 cm 

dz (1 + z) 2 (1 + 10 5 n 0 z) 


(145) 


n 0 being the present mean atomic density of all the matter in the universe. We 
will consider here two types of model universes: 1) a "flat" or Einstein-de 
Sitter model with n 0 ~ 10" 5 cm 3 and 2) an "open" model with n 0 « 10“ 5 z. 
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For the flat model, equation (144) reduces to 


t(E^, , z ) = 3.9 x 1 

v yO max ' 


o 8 f"“ d, »p y (146) 

Jq (1 + z) 1/2 


with E y0 in GeV. 


For z fflaJ£ >> 1, equation (146) can be further simplified to yield 


T <>V z „„)~ 17 x (1 + exp- 


1.12 x 10 “ 

(! + z ma J 2 Ero 


(147) 


A numerical solution found by setting equation (144) for (E , z crit ) = 1, 
which defines the critical redshift where the universe becomes opaque to gamma- 
rays of energy E can be well approximated by the expression 


1 + z crit Of 2,60 x 10 2 E ^ -484 


( 148 ) 


For the open model, we find 


I III O A. 

dz exp 


1.12 x 10 6 

(i + *) 2 K 


( 149 ) 


For z max » 1, equation (149) may be approximated by 


z.,,) 11 - 7 * 102 Ki 1 <1 + *.„> 


exp - 


1.12 x 10 6 

(i + z ) 2 r 


( 150 ) 


Thus, there is no significant difference between the opacities of the open 
and flat model universes. This being the case, we may invert equation (148) to 
obtain an expression forlhe predicted cut-off energy, E c , above which gamma- 


rays originating at a redshift, z max , cannot reach us. This relation is then given 
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E 


2.60 


(151) 


C 


'2.6 Ox 10 2> 

1 + 


and is graphed in Figure 26. 

In the other extreme, v « 1, we find that as we consider higher and higher 
energies, the universe will not become transparent to gamma-rays again until 
we reach an energy E tr , where the optical depth, r (E tr , z max ), again falls to 
unity. The expression for the optical depth when v « 1 in equation (63) is given 
for a flat universe by 


t(E^, z ) ^ 4.4 x 10 13 E 

v yu mflx ' 


p 

r o 

*'0 


dz (1 + z) 


- 3/2 


(152) 


— 8.8 x 10 13 e;J(1 


+ z 


) 1/2 for 


» 1 


and for an open universe by 


T ( E 7 o- 


z ) - 4.4 x io 13 e;' 

max ' yu 



dzCl + z)" 1 - 4.4 x 10 13 E;J ln(l + z max )(153) 


In both cases we find that E tr > 10 13 GeV so that we may safely assume 
that the universe, due to the blackbody radiation field, is essentially opague to 
gamma-rays of all energies greater than E c . 

Gamma-ray absorption by pair production and Compton interactions with 
intergalactic gas at high redshifts has been examined by Rees (1969) and by Arons 
and McCray (1969). Our discussion here essentially follows theirs. The absorp- 
tion cross section for Compton interactions is given in equations (66) and (67); 
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E r (MeV) 

Figure 26. The redshift at which the universe becomes opaque to photons given as a function of 
observed gamma-ray energy. Gamma-rays originating at all redshifts below the curve can 
reach us unattenuated with the energy indicated. The two curves on the left side of the figure 
are for attenuation by Compton scattering with intergalactic electrons having the densities in- 
dicated and for pair production and are based on the calculations of Arons and McCray (1969). 
The right-hand curve results from attenuation of gamma-rays by interactions with the micro- 
wave blackbody radiation and is based on the discussion of Fazio and Stecker (1970). 




that for pair production in equation (69). The cross sections for these processes 


are graphed in Figure 8. The total absorption cross section above 100 MeV energy 
is roughly constant equals 1.8 x 10" 26 cm 2 . For the case of a constant absorption 
cross section, the optical depth as a function of redshifts, z max , is given by equa- 
tions (144) and (145) as 

dz n<z) "(It) 



n crcQ 
c 



dz 


(1 + z) 
1/1 + 'dz 


(154) 


where 





t = n o"cH7 1 = 1.8 x 1 O'* 3 

c c 0 

a result which was first obtained by Gunn and Peterson (1965). 
For an Einstein-de Sitter universe (fi = 1), 

r(z ) = — t [(1 + z ) 3/ 2 - 1] 

v max ' ^ c max 7 

and for a low density universe (Qz « 1) 

r (z ) = — r [( 1 + z ) 2 - 1] 

v max ' 2 c x max / 


(155) 


(156) 


(157) 
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At lower energies, when the cross section is not constant, the calculation 
is more complex. 

Figure 26 shows the critical redshift for absorption of y -radiation plotted 
as a function of observed energy. At lower energies absorption is due to Compton 
interactions with intergalactic matter, in the intermediate range absorption is 
due to pair -production interactions with intergalactic matter (Arons and McCray 
1969 Rees 1969). At the higher energies absorption is due to pair production 
interactions with blackbody photons (Fazio and Stecker 1970). As one can see 
from the figure, there is a natural "window" between ~1 MeV and ~10 GeV 
which defines the optimal energy range for studying cosmological y -rays. 

In the case of pair production interactions, i.e., those which are important 
for > 1 MeV as shown in Figure 26, the photon completely disappears and the 
effect is a true absorption effect. In those cases, particularly above 10 GeV, z c 
can be used as an upper limit on the integral given by equation (105) for evalu- 
ating the background energy spectrum. This immediately implies a steepening 
of the background spectrum above 10 GeV. However, in the case of Compton 
scattering, important for E < 1 MeV, we do not have a pure absorption process 
but rather a process in which a photon of some initial energy E' is replaced by 
one at some lower energy E. Thus, in order to calculate the theoretical back- 
ground spectrum properly one must solve an integrodifferential scattering 
equation (Arons 1971 a, b). We will refer to this equation as the cosmological 
photon transport (CPT) equation. We can write this equation in the form 


88 



M + J_[-EH(z)S]=2<E. «)-*„(E,z) f"’ dE'« sc (E,z»(E|E')dE' (158) 

B 1 BE J E 

where E is the photon energy k ^ and k sc are the photon absorption and scatter- 
ing rates (which are a function of z because the intergalactic gas density is as- 
sumed to scale as (1 + z) 3 because of the expansion of the universe. The script 
quantities for the y -ray intensity and production rate 


<SL (E, z) = (1 + z) -3 I(E, z) 


(159) 


2(E, z) ■ (1 + z)- 3 Q(E. z) 

are quantities co-moving with the expansion, defined so that their redshift-density 
dependence cancels out 8(E) is an upper limit on the scattering integral defined 
by the Compton process and H(z) is the Hubble parameter which, in terms of the 


Hubble H., is given by the relation 


H(z) = H 0 (l + z) /TToi 


(160) 


where H is the ratio of the mean gas density in the universe to the density needed 
to close the universe gravitationally. The term 


B<9. 

Bt 


(1 + z)H(z) 


2z 


(161) 


and the second term in equation (158) expresses the energy loss of the 7 -rays 
because of the expansion redshift. 

The third term in equation (158) is the source function or production function 
of y -radiation produced with energy E at redshift z. The forth term represents 
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the true absorption due to pair -production processes as well as the total scatter- 
ing cross section for Compton scattering as a function of redshift and energy. 
Because all of the scattered photons "disappear" at energy E' but reappear at 
some lower energy E, the integral term in the right-hand side of the equation has 
to be included. This term is equivalent to a source term of photons at lower 
energies. The upper limit on the integral is determined by the kinematics of 
Compton scattering and is given by 

fE/bn c 2 (l - 2E/m c 2 )], E<- m c 2 

S(E) = J (162) 

l 00 . E^Im e c 2 

The CPT equation was solved numerically by Stecker, Morgan and Bredekamp 
(1971) in their study of the spectrum of background-radiation to be expected from 
cosmological matter -anti matter annihilation. Part of their results are shown in 
Figure 25. The curve labeled "no absorption" shows the result with the terms 
involving k ^ and k sc set equal to zero. The curve labeled "absorption without 
transport" shows the result when the integral term involving * sc is set equal to 
zero. This is equivalent to removing all of the scattered photons and not re- 
placing them at lower energies. The solid line shows the full solution. 

Figure 26 shows the effect of increasing the mean gas density in the inter- 
galactic medium and therefore the mean density of Compton scattering electrons. 

It clearly demonstrates the increase in the absorption due to Compton scattering 
at low energies as the mean density increases. 
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7. Interpretation of the Diffuse y-Ray Background Observations 

It must have been evident very early, in an implicit way that the only source 
of matter large enough to give a significant background of isotropic radiation of 
a truly astronomical nature is the universe itself. Therefore, the connection 
with cosmology has been clearly the prime motivation for interest in an isotropic 
background radiation since Morrison's 1958 paper. 

Theoretical discussions of various types may be found in the literature long 
before the first solid observational evidence of the existence of a gamma-ray 
background. 

The mechanisms listed by Morrison (1958) to be of possible significance in 
producing continuum radiation were synchrotron radiation, cosmic-ray electron 
bremsstrahlung, 77°-meson decay (from cosmic ray-nucleon interactions) and 
antimatter annihilation. To these four, one more mechanism, viz., Compton 
interactions between cosmic-ray electrons and starlight (Felten and Morrison, 
1963) was added. It later became apparent when the 3K blackbody radiation was 
discovered that these photons would be orders of magnitude more numerous than 
starlight photons in intergalactic space as targets for cosmic-ray electrons and 
should thus be the prime Compton radiation generators of cosmological interest. 
That the significance of this was readily grasped is obvious from the plethora of 
independent suggestions made immediately after the discovery of the microwave 
background radiation (Felten 1965, Gould 1965, Hoyle 1965, Fazio, Stecker and 
Wright 1966, Felten and Morrison 1966). The relative weakness of intergalactic 
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magnetic fields, evidenced by data on the non-thermal radio background, elimi- 
nated synchrotron radiation as a prime contender in generating the diffuse gamma- 
ray background so that four mechanisms were left 

(a) electron bremsstrahlung 

(b) electron-photon interactions (Compton effect) 

(c) cosmic-ray produced 77 “-decay 

(d) decay of tt “-mesons from matter -antimatter annihilations 

Evidence for a diffuse background above 50 MeV was reported by Kraushaar 
and Clark (1962) from measurements on Explorer 11. The interpretation of 
Felten and Morrison (1963) that both the Ranger 3 and Explorer 11 results could 
be fitted reasonably well by a single power law of the type expected from Compton 
interactions seemed logical despite a possible flattening above 1 MeV reported 
by Arnold, et al. (1962) and it was expected that data in the two energy decades 
between 1 and 100 MeV would exhibit nothing more exciting than a smooth power 
law spectrum as extrapolated from the sub-MeV ("X-ray”) energies. 

In the late 60's the author, after having made detailed thesis calculations of 
gamma-ray spectra from cosmic-ray producted secondary particles and from 
proton-antiproton annihilation (Stecker 1967), became interested in the effects 
cosmology might have on such spectra and on the implications of these effects 
for cosmology itself. Cosmic-ray 77 ° -decay was suspected to play a major role 
in generating galactic gamma-rays (Pollack and Fazio 1963) and it remained a 
viable possibility the extragalactic background flux above 100 MeV. But if such 
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interactions are occurring in intergalactic space now, why not in the distant 
past when gas and cosmic-ray densities were higher (in an expanding universe)? 

If so, large fluxes of extragalactic cosmic-rays (comparable to galactic fluxes) 
need not exist now to explain the 100 MeV background (Stecker 1968, 1969a). 

Also, the spectrum would be redshifted and would be softer than the galactic 
spectrum (Stecker 1969b). Similar ideas were being independently worked on 
by Ginzburg (1968) in the context of the Lamaitre cosmological model and by 
Rozental and Shukalov (1969) for the standard expanding universe model. In 
these models, various cosmological effects come into play to distort the spec- 
trum from t r ° decay and redshift its characteristic peak from an energy of 
m^c 2 ^ - 70 MeV to lower energies (see extensive discussions in the mono- 
graphs by Stecker (1971a) and Ozernoi, Prilutsky and Rozental (1973)). The 
result was the prediction of a possible enhancement in the gamma-ray background 
spectrum between 1 and 100 MeV deviating from the simple power law extra- 
polation of the X-ray background. 

Figure 29 shows a two component model normalized for a best-fit to the 
observations involving the production of intergalactic gamma-rays from cosmic- 
ray interactions with intergalactic gas producing 77 “-mesons out to a maximum 
redshift of 100 (Stecker 1969b, c, 1971b). Three problems arise with this ex- 
planation: 1) even with a relatively steep assumed cosmic-ray spectrum 
(~ E ~ 2- 1 ) the bulge in the theoretical spectrum may be too large to fit the ob- 
servations, although this discrepancy may not be too serious considering 
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Figure 27. The cosmological 7 -ray spectrum from matter- 
antimatter annihilation calculated by solving the CPT equa- 
tion numerically for Q =1. The solid line represents the 
complete solution. The other curves represent the effect of 
neglecting the absorption and scattering (transport) terms 
in equation (158). 
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ARBITRARY NORMALIZATION 



E r (MeV) 

Figure 28. The effect of absorption of y -rays at high redshifts 
by the protogalactic gas. 
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I(Ey) (cm' s.$r.MeV)‘ 



Figure 29. Comparison of the observed background 
with a two-component model involving the produc- 
tion and decay of neutral pions produced in inter- 
galactic cosmic -ray interactions at red-shifts up to 
100 . 


96 



observational uncertainties, 2) large amounts of energy are needed in cosmic- 
rays at high-red-shifts, requiring the existence of strong primordial cosmic-ray 
sources (protars) which are either pregalactic or protogalactic. Such objects 
may have been the "spinars" suggested by Morrison (1969) if spinars existed 
at such redshifts of about 70-100 (Stecker 1971b). If it is considered that each 
spinar produces approximately 10 62 ergs over a time scale of 10 7 -10 8 years 
(Morrison 1969), a time comparable to the Hubble time at these redshifts, then 
at most 20 percent of the presently observed galaxies are needed to have arisen 
from this early spinar state in order to provide the cosmic-ray energy needed 
to account for the diffuse y-radiation above 1 MeV. At a redshift of about 70, 
the free-fall time for forming spinars from gas clouds is comparable to the 
Hubble time. This may provide a natural upper limit to the redshift, z,^ , for 
primordial cosmic-ray production in the spinar model. (It should, however, be 
noted that such spinars may arise in other ways (see Stecker 1971b) and that 
they may now be a class of moribund objects unrelated to galaxies as we see 
them now). 3) The third problem with this hypothesis is that the maximum red- 
shift for cosmic ray production z max is a free parameter chosen to fit the 
observations. 

A related hypothesis examined by Stecker et al. (1971) and Stecker and Puget 
(1972) is that the y -ray background is from redshifted 7r°-decay y-rays but that 
77 °-mesons are the result of nucleon-antinucleon annihilation at an early epoch 
in the history of the universe. The annihilation hypothesis does not suffer from 
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the above mentioned problems of the cosmic-ray protar hypothesis. The param- 
eter z max does not enter into the theory; annihilations occur at all redshifts and 
the 1 MeV-flattening is an absorption effect as discussed earlier. The transport 
equation (158) was solved to determine the exact form of the spectrum. Energy 
considerations do not present a problem. Another advantage of the theory is that 
it arises as a natural effect in a cosmology such as that suggested by Omnes 
(1972 and references therein) Figure 30 shows a detailed comparison of the 
annihilation hypothesis spectrum with present observations. The dashed line 
shows the effect of adding in an additional component which is a power -law 
extrapolation of the x-ray background. The two component model shown provides 
an excellent fit to the observational data. 

Early observations of gamma-radiation by Kraushaar and Clark (1962) had 
clearly indicated that if antimatter exists in the universe in large amounts, it 
must clearly be separated from matter so that the average annihilation rate is 
quite small. In 1969, Omnes suggested a baryon symmetric (equal amounts of 
matter and antimatter) cosmology based on a possible phase transition effect 
which could separate matter from antimatter at an early stage in the big-bang 
corresponding to nuclear density for the cosmic plasma. The phase transition 
effect was also studied by Aldrovandi and Caser (1972) and Cisneros (1973). 
Further work by Omnes (1972 and references therein) showed that the separate 
domains of matter and antimatter could grow to contain masses of the size of 
galaxies by the recombination epoch. This result has recently been refined by 
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Figure 30. A comparison of the data given in Figure 25 
with the annihilation model discussed by Stecker, 
Morgan and Bredekamp (1971) and Stecker and Puget 
(1972). 
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Aldrovandi et al. (1973). It was to be expected that boundary-region annihilations 
in this picture would also produce redshifted 7r°-decay radiation and absorption 
effects would cut off the resultant flux below 1 MeV. Therefore, Stecker, Morgan 
and Bredekamp (1971) were motivated to make a detailed calculation of the re- 
sultant diffuse background spectrum to be expected, and the results agreed fairly 
well with the observations then available. The encouraging enhancement in the 
1 to 100 MeV range is partially due to the existence of a "gamma-ray window" 
in this energy range as shown in Figure 26. The results were encouraging enough 
to examine further the evolution of the Omnes cosmology for redshifts less than 
10 3 (Stecker and Puget 1972). This study had several exciting implications 

(a) separate regions containing masses the size of galaxy clusters could 
be obtained. 

(b) turbulence produced by annihilation pressure could provide enough 
energy to trigger galaxy formation, 

(c) estimates obtained placed the galaxy formation stage at redshifts of 
the order of 60. 

(d) mean densities and angular momenta of galaxies could be estimated in 
this picture consistent with observation and related to the annihilation 
rates calculated by the model and implied by the observations. 

The general scheme of the galaxy formation model is shown in Figure 31. 

The observational implications of the model are outlined in Figure 32. 
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GENERAL SCHEME OF MODEL 


• matter and antimatter exist in equal amounts in separate regions. 

• MIXING OCCURS ALONG BOUNDARY REGIONS OF THICKNESS~A a . 

• RESULTING RAPID HEATING OF PLASMA WITHIN A DISTANCE ~ Ax OF BOUNDARY BY 

ANNIHILATION PRODUCTS PRODUCES EXPANSION AWAY FROM BOUNDARY. 

• RESULTING EXPANSION OF PLASMA INDUCES HIGH-VELOCITY GAS MOTIONS. 

• DURING THE NEUTRALIZATION ERA (POSSIBLY EVEN SOMEWHAT BEFORE) GAS MOTIONS BECOME 

TURBULENT (~500<Z <~3000). 

• WHEN GAS GOES FROM PLASMA TO ATOMIC STATE <~400 < Z N <~600) : 

A) DUE TO DECOUPLING OF MATTER FROM RADIATION FIELD, VISCOSITY DROPS BY 
ALMOST 8 ORDERS OF MAGNITUDE AND SOUND VELOCITY DROPS BY FOUR ORDERS 
OF MAGNITUDE. 

B) THIS CAUSES TURBULENCE TO BECOME SUPERSONIC AND TO EXTEND THE EDDY 
SPECTRUM DOWN TO A SCALE OF THE ORDER OF 10' 3 pc. 

C) THE SUPERSONIC TURBULENCE INDUCES DENSITY FLUCTUATIONS Ap/p~10VER 
THE WHOLE RANGE OF EDDY SCALES. 

• ATAREDSHIFTOF ^-Z N OR ABOUT 60, VIRIAL THEOREM BECOMES SATISFIED FOR BINDING OF 

GAS CLOUDS INTO PROTOGALAXIES DUETO DENSITY FLUCTUATIONS INDUCED BY 
SUPERSONIC TURBULENCE. 


Figure 31. Outline of the galaxy formation theory of Stecker and Puget (1972). 
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Figure 32. Observational implications of baryon symmetric cosmology. 




















Several other models of isotropic y-ray production have been put forward. 
One suggestion is that the whole spectrum in the 10 - 3 -10 2 MeV range is due to 
Compton interactions of intergalactic electrons with the universal blackbody 
radiation (Felton 1965, Gould 1965, Hoyle 1965, Fazio, Stecker and Wright 1966, 
Felten and Morrison 1966). In its most recent version Brecher and Morrison 
(1969) have attempted to explain the observed spectral features using the Compton 
hypothesis, however, Cowsik and Kobetich have done a more detailed calculation 
indicating that the Compton mechanism generates a smooth featureless power- 
law spectrum. 

Figure 32 shows the energy spectrum J(E y ) = E y I(E y ) of the background 
radiation between 10“ 3 and 200 MeV as based on the review paper of Schwartz 
and Gursky (1973) and the data shown in Figure 23. An extrapolation of the data 
between 30 keV and 1 MeV is shown by the dashed line. A strong deviation from 
a power -law spectrum is indicated. 

A critical test between the cosmic ray proton and annihilation hypotheses 
lies in a study of the energy spectrum. Figure 33 shows the present range of 
the data, indicated by the shaded region, along with the extrapolated power-law 
spectrum (X) the annihilation spectrum (A) and the high energy form of the 
spectrum predicted for redshifted cosmic-ray 7 ?° -decay gamma-rays (CR). 

The annihilation spectrum should exhibit a sharp cutoff slightly below 1 GeV 
because the energy of the gamma-rays is limited by the rest-energy available 
to them from baryon-antibaryon annihilations. A detailed discussion of this 
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E 7 (MeV) 

Figure 33. Energy flux spectrum of the x-ray and gamma-ray background 
based on Schwartz and Gursky (1973) and the data given in Figure 25. The 
straight diagonal line indicates an extrapolation of the 30 keV to 1 MeV 
spectrum. 
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may be found in Stecker (1971a). The cosmic-ray produced spectrum, on the 
contrary, can continue up to higher energies with a steepening induced around 
10 GeV by pair -production losses through interactions with the microwave back- 
ground (Fazio and Stecker 1970). This steepening should amount to an increase 
of 0.5 in the spectral index for a closed Einstein-de Sitter universe and an in- 
crease of 0.75 in the spectral index for a low-density open universe (Stecker 
1971a). It should be kept in mind that the cutoff in the annihilation spectrum 
may be somewhat obscured by the presence of other background radiations having 

relatively lower intensities below 200 GeV. 

Various other mechanisms which have been discussed in connection with the 
interpretation of the origin of the y-ray background may be found in a review last 
year by the author (Stecker 1973c). 

Another recent hypothesis for the origin of the background radiation which 
also involves primordial cosmic rays is that of Strong, et al. (1973). This 
hypothesis is based on the model of Hillas (1967) that the observed cosmic rays, 
at least those above ~10 6 GeV, are universal and primordial. Hillas showed 
that if the observed cosmic rays originated at a redshift z max - 15, energy losses 
that these cosmic rays would undergo in pair -production interactions with the black- 

body photons of the type 

y + p— p+e + +e" (163) 

would steepen the spectrum at the presently observed energy of ~ 3 x 10 6 GeV. 
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At a redshift of ~15, when these interactions took place, the blackbody photons 
had an energy of ~ 10" 2 eV and the cosmic rays had an energy of ~5 x 10 7 GeV. 
The electron-positron pairs produced would then produce a shower of lower 
energy 7 -rays through a repeated two -stage interaction process with the black- 
body radiation field involving Compton scattering followed by pair production, 
i.e., 


+ ^blackbody e + ^high energy 

followed by (164) 

^blackbody ^Hi gh energy e 4 e 

The resultant spectrum based on a complex model is shown in Figure 35. 

Both the Hillas model and the hypothesis of Strong, et al. require an exceedingly 
low intergalactic gas density, n 0 < 10" 9 cm -3 , in order that the y -rays produced 
by the cosmic ray model of Stecker (1969b) do not exceed the observed background 
level. 

If the intergalactic gas density should turn out to be of the order of 10“ 9 cm -3 
which is only about 1% of the mean matter density in galaxies, we would have to 
radically alter our ideas about galaxy formation in order to account for 99% of 
the mass in the universe to be bound into galaxies. 

Present observations give a background flux of y radiation above 100 MeV 
of about 


J 


qL 

4tt 


3x 1 0 -5 cm -2 s _1 Sr -1 


(165) 
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Figure 34. Predicted energy flux spectra from the annihilation model 
(A) and cosmic ray (protar) model (CR) as discussed in the text. Also 
shown is the scatter area covered by the observational data (shaded) 
and the extrapolated x-ray background spectrum (X). The two curves 
shown for the CR spectrum above 7 GeV are for closed (Einstein-de 
Sitter) and open universes as discussed in the text (Stecker 1971a). 
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Figure 35. The 7 -ray background calculated from the Hillas 
model by Strong, et al. (1973) (SW 2 ). Also shown are the 
data from Figure 25 for comparison. 
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where the path length 


L - - 10 28 cm 

H-o 


(166) 


and the production rate 


q - 1.3 x 10" 25 


(167) 


2? io -30 n<f 


where <f is the ratio of intergalactic to galactic cosmic ray intensity if part of 
the background is produced by noncosmological (contemporary z = 0) cosmic 
rays. If we call the fraction of the background produced by contemporary cosmic 
rays f, then from equations (165) to (167) 

3 x 10 -s f 2? 10' 3 Of ( 168 ) 


because the background spectrum is observed to be quite steep, conservatively 
f < 1/3 so that the limit on the product of intergalactic gas density and cosmic 
ray intensity becomes 

Q£< 10' 2 ( 169 > 

This excludes the possibility of both universal cosmic rays (f = 1) and a closed 
universe (fl = 1). 

More stringent upper limits can be placed on primordial cosmic rays 
(Stecker 1971b) and the antimatter annihilation rate (Stecker, et al. 1971), Stecker 
and Puget 1972). Thus it would seem that y -ray astronomy at present can with 
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certainty only put limits on our speculations about the origin of cosmic rays. 
But on the other hand we can say that speculations about the origin of cosmic 
7 -rays have broadened our conceptions about high-energy astrophysics and 
cosmology. 
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